Phase separation in a polarized Fermi gas at zero temperature 
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We investigate the phase diagram of asymmetric two-component Fermi gases at zero temperature 
as a function of polarization and interaction strength. The equations of state of the uniform su- 
perfluid and normal phase are determined using quantum Monte Carlo simulations. We find three 
different mixed states, where the superfluid and the normal phase coexist in equilibrium, corre- 
sponding to phase separation between: (a) the polarized superfluid and the fully polarized normal 
gas, (b) the polarized superfluid and the partially polarized normal gas and (c) the unpolarized 
superfluid and the partially polarized normal gas. 

PACS numbers: 



The study of a two-component Fermi gas with imbal- 
anced populations is an active area of research in the 
field of ultracold atoms Recent experiments, carried 
out on harmonically trapped configurations, investigate 
superfluidity and Bose-Einstein condensation (BEC) of 
fcrmionic pairs in these systems by varying the strength 
of interactions, the temperature of the gas and the de- 
gree of polarization Evidences of phase separation 
between a superfluid core and a normal external shell are 
reported for systems close to the unitary limit 0] and on 
both sides of the Feshbach resonance Q ■ On the theoret- 
ical side, the phase diagram of a polarized Fermi gas with 
tunable interactions is the subject of a number of studies 
both at zero [B|, Q and at finite temperature 0] • These 
studies, which refer to uniform systems and are based on 
a mean-field approach, predict the existence of a normal 
phase for large enough polarization on the BCS side of 
the crossover and of a polarized superfluid phase on the 
BEC side, separated by a region where the two phases 
coexist in equilibrium. However, if interactions are not 
weak, the mean-field theory fails to describe correctly the 
nature of the phase separated state and to provide a re- 
liable estimate of the critical value of polarization where 
phase separation occurs. The key ingredient, that is miss- 
ing in the mean-field description, is the proper account 
of interaction effects in the normal phase Q. 

In this Letter we carry out a quantitative study of the 
phase diagram at zero temperature along the BCS-BEC 
crossover using fixed-node diffusion Monte Carlo (FN- 
DMC) simulations, which have proven very accurate in 
the investigation of strongly correlated ultracold Fermi 
gases [l[ . We determine the equation of state of the ho- 
mogeneous normal and superfluid phase as a function of 
interaction stregth and population imbalance. From the 
phase equilibrium conditions we identify three different 
phase separated states corresponding to: (a) the polar- 
ized superfluid coexisting with the fully polarized nor- 
mal gas, (b) the polarized superfluid coexisting with the 
partially polarized normal gas and (c) the unpolarized 
superfluid coexisting with the partially polarized normal 
gas. The state (a) and (c) are respectively separated from 
the homogeneous superfluid and normal phase by a first 



order phase transition, while state (a) and state (b) as 
well as state (b) and state (c) are separated by second 
order phase transitions. 

We consider a uniform system in a volume V with a 
total number of particles N — + iVj,. The number A/j 
of spin-up particles is kept fixed, corresponding to the 
average density — iVj /V. The interaction strength 
is parametrized by the inverse product l/fc^-j-a of the s- 
wave scattering length a and the Fermi wavevector of the 
spin-up particles fe^-j = {fy^n^) 1 / 3, . The number iVj of 
spin-down particles is instead a variable of the system de- 
termining the polarization parameter P = +N l > that 
is assumed positive within the bounds < P < 1. The 
relevant energy scale is fixed by the Fermi energy of the 
spin-up particles Epi = h 2 kp^ /2m, where m is the par- 
ticle mass of both spin components. 

First, we discuss the equation of state of the four ho- 
mogeneous phases considered in the present study: un- 
polarized and polarized superfluid and fully and partially 
polarized normal gas. Then, we analyze the equilibrium 
conditions between these phases to map out the phase di- 
agram along the BCS-BEC crossover. We do not consider 
more exotic superfluid phases such as the Fulde-Ferrel 
and Larkin-Ovchinnikov (FFLO) state which is expected 
to occur on the BCS side of the resonance for small po- 
larizations Q. 

i) Fully polarized normal gas (Npp). 

In this phase Nt = 0. Since p-wave collisions can be 
neglected, the gas is well described by the non-interacting 
model with the equation of state 



(1) 



ii) Unpolarized superfluid gas (SFo) 



In this phase Nt 



N/2. The corresponding 



ground-state energy can be cast in the following form 



N 



SFo 



e 6 6(l/fc FT a)- 



-N t E 



2G(l/fc FT a) . (2) 



The first term, proportional to the number of pairs N/2, 
corresponds to the contribution of molecules with bind- 
ing energy For a zero-range interatomic potential this 
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FIG. 1: Binding energy A of a single spin-down impurity 
in a Fermi sea of spin-up particles. The dashed line is the 
molecular binding energy et for our short-range square well 
potential (with 2ti^Rq = 10~ 6 as in Ref. [HI). In the inset 
we show the equation of state of the unpolarized superfiuid 
SF . The solid lines correspond to best fits to the FN-DMC 
results. 
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FIG. 2: (color online). Equation of state of the normal par- 
tially polarized phase Npp as a function of the concentration 
x for different values of the interaction strength. The solid 
lines correspond to best fits of the energy functional © with 
the values of A and m* obtained from the single-impurity 
calculations. 



energy takes the familiar expression = —h 2 /ma 2 . The 
two-body term is a convenient parametrization of the 
equation of state only on the BEC side of the crossover 
(a > 0) where dimers are formed in vacuum, as en- 
tailed in Eq. @ by the Heaviside function: Q(x) = 1 
if x > and zero otherwise. The dimensionless function 
of the interaction strength G(l/kF]Ci) contains instead 
the many-body contributions to the equation of state. 
This function has been calculated in Ref. [lfj using the 
FN-DMC method and the results are shown in the inset 
of Fig. [TJ In particular, in the BEC regime ^S> 1) 

the function G takes the form G = 5fc ft add /187r[l + 
128(^ T a d(i /7r) 3 / 2 /15\/6] with a dd = 0.60a [Hf and corre- 
sponds to the mean-field and first beyond mean-field con- 
tributions to the equation of state of composite bosons 
with mass 2m and density nt interacting with a dimer- 
dimer scattering length add- 

Hi) Partially polarized normal gas (Npp). 

This phase is characterized by the concentration x = 
iVl/iVf °f the minority spin-down particles. At small con- 
centrations (i « 1) the dependence on x of the ground- 
state energy can be written in the form of the Landau- 
Pomcranchuk Hamiltonian of weakly interacting quasi- 
particles [12j 

£ Npp = ^N,E F1 ( 1 - Ax + — x 5/3 + Fx 2 ) , (3) 
5 V to* / 

where A, m*/m and F depend on l/kp\a. The quan- 
tity A is the binding energy of a spin-down quasiparti- 
cle in the Fermi sea of spin-up particles and to* is its 
effective mass. The term F accounts instead for the cou- 
pling between quasiparticles. As already pointed out, 
the mean-field approach completely neglects interactions 



in this phase resulting in an energy functional given by 
Eq. ([3]) with A = F = and to*/to = 1. We determine 
A and to* as a function of the interaction strength from 
quantum Monte Carlo simulations. The binding energy 
is obtained from the ground-state energy of the system 
with one spin-down impurity in a Fermi sea of spin-up 
particles and the effective mass from the curvature of the 
excitation energy if the impurity carries a small momen- 
tum The results for A are shown in Fig. [TJ They are 
in excellent agreement with the values recently obtained 
using exact diagrammatic Monte Carlo methods [l]| and, 
quite remarkably, also with the results of a simple varia- 
tional approach based on a single particle-hole wavefunc- 
tion [lj] . The values we obtain for the effective mass m* 
are instead slightly smaller than the diagrammatic Monte 
Carlo results of Ref. [13| . This might be due to a non 
optimal choice of the nodal surface of the excited state 
at finite momentum, the FN-DMC method provides in 
fact only an upper bound for the energy unless the nodes 
of the many-body wavefunction are known exactly, and 
to finite size effects in the analysis of the low-momentum 
spectrum [15| . We also perform FN-DMC calculations at 
finite concentration x for various values of the interaction 
parameter l/fc^a using the Jastrow-Slater wavefunction 
described in Ref. [l2[. The results are presented in Fig. [5] 
By fitting the functional form ([3]) to these results with A 
and to* obtained from the single- impurity calculations, 
we determine the interaction parameter F and its de- 
pendence as a function of 1/kp^a. At unitarity we find: 
A(l/k F] a = 0) = 0.99(1), TO*(l/fc FT a = 0)/to = 1.09(2) 
and F(l/kF\CL_= 0) = 0.14, in agreement with the find- 
ings of Ref. [12j . We notice that the Npp phase reduces 
to the Npp one if x = 0. 

iv) Polarized superfiuid gas (SFp). 
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FIG. 3: (color online). Equation of state of the superfluid 
polarized phase SFp as a function of the concentration y for 
different values of the interaction strength. The solid lines 
correspond to the energy functional Q. 



This phase is characterized by a number Np — iVj_ of 
pairs and a number Na = N-\ — N± of unpaired atoms, 
such that N = 2Np + Na- We denote the concentration 
of the minority atoms by y = N± /N-\ . In the deep BEC 
regime the SFp phase corresponds to a miscible mixture 
of Np bosons and Na fermions 

EE 

The interaction 

between bosons and fermions is repulsive and is fixed by 
the atom-dimer scattering length a a d — 1.18a |18j |. In 
this regime we write the equation of state in the form 



Esf (Np 

— iV| Ep^ 
5 



(4) 



(i-y) 5/3 + 



5k F ia ad 
3tt 



2/(1 - y) 



where E SFg {N P ) = 3/5iV T £; FT2 ; 5 / 3 2G(l/fc FT ay 1 / 3 ) + 
e^Np is the energy of the 2Np paired atoms and the 
other terms in Eq.((3]) correspond to the kinetic energy 
of the unpaired atoms and to the interaction energy be- 
tween atoms and dimers treated at the mean-field level. 
We carry out FN-DMC simulations of the SFp phase for 
various values of the interaction strength. The nodal 
surface is modeled using a BCS plus unpaired particles 
wavefunction written in the form of a determinant as in 
Ref. [l9j]. The results are shown in Fig. [3] together with 
the energy functional (TJJ. The agreement is remarkable 
down to quite small values of the interaction parameter 
l/kF^a > 0.5[20(. Furthermore, we notice that the SFp 
phase reduces to the SF one if y = 1. 

We are now in a position to study the coexistence 
between the superfluid and normal phases introduced 
above. One requires the equilibrium of pressures be- 
tween the superfluid ps — —dEs/dVs and the normal 
Pn = —dEpf/dVjsr state and the equilibrium of chemical 
potentials. In the normal phase there are two chemi- 
cal potentials for the Npp state: MjvtU) = dEx/dN^m, 
which reduce to only /in^ for the Npp state. Similarly 
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FIG. 4: (color online). Phase diagram as a function of po- 
larization and interaction strength. In terms of the Fermi 
wavevector fcf = (3n 2 n) 1 ^ 3 fixed by the total density n = 
N/V one has: 1/kpa — 1/kFia at P = and I/Zcfh = 
2 li/3 /fcFT a at P = 1. On the BCS side of the resonance our 
determination of the critical polarization is not reliable. For 
— \/kF]CL > 1 we obtain P c using the BCS theory (see text). 



in the SFp phase one can vary both the number of pairs 
Msp = dEs/dNp and the number of unpaired atoms 
Msa — dE$ /8Na, while in the SFo phase only the chem- 
ical potential of pairs /isp with Np = N/2 is relevant. 

a) Phase separation between SFp and Npp. The equi- 
librium conditions are: ps — Pn and pL$A — ^wt- F° r 
a given concentration y of the spin-down atoms in the 
SFp phase, the two conditions determine the values of 
the densities of the spin-up component in the two co- 
existing phases. The ratio P c = gives the critical 
polarization above which the system begins nucleating 
the normal phase to accommodate the excess polariza- 
tion. By increasing P above P c the equilibrium densities 
of the two phases, as well as the critical concentration 
y of the SFp phase, do not change, instead, the vol- 
ume fraction Vn /V of the normal phase increases and 
eventually becomes the entire volume for P = 1. The 
critical polarization line, corresponding to a first order 
phase transition, is shown in the phase diagram of Fig. 2] 
At P = 1 this line terminates at the tricritical point 
l/kp^a = 1.7 |21(. For larger values of l/kp^a the ho- 
mogenous SFp phase exists up to P = 1 and the super- 
fluid to normal transition becomes second order. For a 
given concentration y the SFp-Nfp state is stable pro- 
vided pisp < f+N]- + f-iNl = Epi(l — 3^4/5), i.e. until 
the process in which pairs break and spin-down parti- 
cles start to populate the normal phase remains energet- 
ically unfavourable. The instability line, corresponding 
to a; = 0, marks a second order phase transition where 
the fully polarized evolves continously into the partially 
polarized normal phase. At P — 1 this line terminates at 
the point 1/kp^a = 0.73, for smaller values of l/fc^a a 
superfluid can not exist up to P — 1 . For small polariza- 
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tions the second order transition line terminates at the 
point corresponding to P c — 0.015 and 1/fcpfa = 0.61. 

b) Phase separation between SFp and Npp. In this case 
one has to fulfill three equilibrium conditions: ps = Pn, 
Usa = A*iVT an d fJ-SP = Mjvt + M-ZVl- It is worth point- 
ing out that the SFp-Npp phase separated state does not 
exist within the mean-field description, where either the 
normal state is fully polarized or the superfluid state is 
unpolarizcd. By approaching the unitary limit the SFp- 
Npp state becomes unstable because the polarization of 
the superfluid is energetically unfavourable and the SFp 
phase evolves continously into the SFo phase. The insta- 
bility line, corresponding to y = 1, marks another second 
order phase transition. This line terminates at the fol- 
lowing points: P = 0.94 and = 0.63 and P = 
and l/fc^a — 0.53. In particular, the point at P = 
corresponds to the smallest value of l/fc^-pa below which 
the superfluid phase can not be polarized [22j]. The two 
second order transition lines at x — and y = 1 are 
indeed very close and, given the uncertainty in the deter- 
mination of the equation of state of the various phases, 
we can not exclude that they might coincide, correspond- 
ing to a single second order phase transition where the 
superfluid polarizes and the normal phase becomes fully 
polarized, or that they might have a reversed order pro- 
ducing a small region of SFo-Npp mixed phase instead 
of the SFp-Npp one. It is important to stress, however, 
that in the relevant region we can provide a reliable de- 
scription of the equation of state of the superfluid (see 
Fig. [3]) and of the normal phase, where for small concen- 
tration x only the term containing the binding energy A 
is important. 

c) Phase separation between SFo and Npp. One has to 
fulfill two conditions: ps = Pn and ^sp = Mnt + A*JVJ.- 
The SFo-Npp state is stable provided /ijvt ^ fJ-SAi cor- 
responding to the instability against polarization of the 
superfluid. The instabilty line coincides with the second 
order phase transition for y = 1 obtained above. The 
critical polarization line P c = marks instead the first 
order phase transition from the normal to the superfluid 
gas. At unitarity we find P c = 0.39 [12j in contrast with 
the value P c = 0.93 predicted by mean-field theory 0. 
Notice that the Landau-Pomeranchuk energy functional 
([3]) does not provide a valid description of the Npp phase 
if the concentration x becomes large. For this reason, on 
the BCS side of the resonance, our results for P c are lim- 
ited to the region close to the unitary limit. In the deep 
BCS regime one can calculate the critical polarization by 
using the following energy functionals for the Npp and 
SFo phase respectively 
20A)Ft|a|a;/97r) and E SF 
where A» a „ 



(2/e) 7 / 3 E F ^ 



HEl: E Npp = 37V T £ FT /5(1- 
#N PP - 3iV T A^ ap /4i; FT , 
■/2k Fr \a\ jg tlie superfluid 

gap. From the equilibrium conditions one obtains P c = 
3/V8(2/e) 7/3 e" 7r / 2fcF Tl Q l, holding to leading order in P. 

In conclusion, we have investigated the phase diagram 
of a Fermi gas at T — as a function of polarization 



and interaction strength. This analysis, carried out for 
uniform gases, is relevant also for systems in harmonic 
traps that can be studied by means of the local density 
approximation. 
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